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Abstract: We investigate a dynamically generated Mott gap from holographic 
fermions in asymptotically geometries with hyperscaling violation by employing a bulk 
dipole coupling for the fermion field. We find that when the coupling strength in- 
creases the spectral function first transfers to the negative frequency region but soon 
redistributes to the positive region. A stable gap and two bands emerges for all momen- 
tum when the coupling strength beyonds a critical value. Generally, The upper band 
on the positive frequency axis is much sharper than the lower band on the negative 
side. When the diploe coupling increases further, the gap becomes larger and the up 
band still keeps sharp while the lower band disperses and widens, concentrating on the 
small momentum region. We also find that the bands will be smoothed out gradually 
with the increasing of hyperscaling violation. 
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1. Introduction 

Recent years, AdS/CFT correspondence has been widely used to study condensed- 
matter theory (AdS/CMT). The strongly coupled conformal theory in the boundary 
is mapped to weakly coupled gravity theory in the bulk. With this great advantage, 
people have successfully constructed holographic models of Fermi and non- Fermi liquids 
in kinds of geometries |L], |2|, ^ f|, ||, [| [7|, 0, [| and can analytically investigate the liquids 
properties, showing the dispersion relation and the width of the quasi-particle like 
excitation. 

Since condensed-matter systems are usually described by non-relativistic field theo- 
ries, to search more proper gravity duals people have further generalized the correspon- 
dence to non-relativistic holography, typically with anisotropic scaling behaviors for 
temporal and spatial coordinates^, [11], [12|, [13| i.e. Lifshitz-like geometry with dynami- 
cal exponent. For realistic systems, importantly, another exponent, called hyperscaling 
violation will emerge and play a crucial role in low energy physics. It is certainly 
necessary to extend holography to this non-trivial case and it indeed comes true by 
employing the standard Einstein-Maxwell-dilation action in the bulk[^, [TJ|, |T5| [HJ |TT| . 
The metric behaves like: 

ds l = -^ + r *n dr 2 + d A (1.1) 
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where i = 1,2, is space index, m and n are related to dynamical exponent z 
and hyperscaling violation exponent 9 by 



m+n+1 n+1 , . 

n + 2 n + 2 

Note when n = —1, the metric reduces to the pure Lifshitz spacetime and n = —2 
corresponds to a class of spacetime conformally related to AdS2 x Rd with the locally 



critical limit z — > oo, 6 — > — oo, while z/6 fixed to be a constant ||18||. The metric 
transforms as 



t ->• A**, ->• Ax,, r -> A ( ^ e)/d r, ds -)• A e/d ds (1.3) 

Clearly, the metric is not scale invariant, unlike the pure Lifshitz case. The dual 
boundary theory exhibits this peculiar behavior generally just below some non-trivial 
dimensional scale. But we won't consider this complication, simply assuming the metric 
is asymptotically geometries with hyperscaling vilation for our purpose. 

We have found remarkable influence of hyscaling violation on the dynamical gap 
by introducing a magnetic dipole coupling for bulk fermions which was first proposed 
in [TJI |20[ and further studied in [^T], [22|, |23|, . Similar to the published work, a gap 
in the spectral function was opened when the dipole coupling strength p exceeds some 
critical value, very like a Mott insulator and further widened with p increasing. The 
coupling strength p plays a role similar to the dimensionless interaction strength U/t 
in the Hubbard model of fermions. The novel feature we find is that two bands exist in 
the spectral function, a upper band on the positive frequency axis and a lower band on 
the negative side, respectively, and behave qualitatively different with the increasing of 
the interaction p and hyperscaling violation 9. 

This paper is organized as follows: In section 2, we briefly review the effective 
gravity model i.e. Einstein-Maxwell-dilation theory for geometries with hyperscaling 
vilation. In section 3, we study the bulk fermions with a dipole interaction, deriving 
the equations of motion for the retarded correlator. In section 4, we numerically solve 
the equations of motion under proper boundary conditions and extract the main results 
of the emergence of the gap. Finally, we present a conclusion in section 5. 

2. Effective Gravity Model 

The standard Einstein-Maxwell-dilaton (EMD) action reads 

S = J d d+2 xV^ [R - 2(<90) 2 - V(<f>) - yZ(0)F 2 - ^H 2 ] (2.1) 
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where the AdS radius has been set to 1. The solutions with hyperscaling violation 
are listed in the following: 

rfs 2 = _ r -2m h ^ dt 2 + r 2n h -l^ dr 2 + ^ ^ ^ = 1 _ ^r_y ^ 

r T'h 

prt = For (m-n+d) z -l^ ? jjrt = ff^m-n+d) ^.3) 



= /c o logr, k = \l^(m-n-2) (2.4) 

\/(0) = -\/ o e- W , K, = 5(m + rf-l), /3 = 2(n ; + 1} (2.5) 

Z-W = Z e-°* + Zl ,a = 2( " + / +1) • % = • * = -f§ (2.6) 

where 5 = m+n+d+1, r\ is the location of the horizon, F Q , H Q are constants which 
are proportional to the conserved charges carried by the black brane. The Hawking 
temperature and the entropy density of the black brane are give by 

T = - 1 = 1 (9 7) 

4 7 r r ( m +"+ 1 ) ' S %K 2 r d h 1 ' ' 

h n 

In the zero temperature limit (r^ — > oo), the entropy density approaches to zero, 
qualitatively differing from the black holes in RN background and more important for 
realistic systems with degenerate ground states. 

In order to admit a stable theory, the dilaton solution is required to be real, leading 
to m > n + 2 or equivalently z > 1 + 9/d, 9 < d. Moreover, in the boundary limit 
r — > 0, the field strength F^ u diverges such that the dual chemical potential cannot be 
well defined. Therefore we introduce another gauge field H = dB to obtain a proper 
definition for finite density 

r (d— m+n+1) 

*(0 = Ml - dt (2.8) 

r h 

where \x is the chemical potential. The constrained condition satisfying B and H 
finite in the UV limit is 

2<m-n<d, d>3 (2.9) 
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The divergent behavior of the field F^ v certainly needs to be treated properly in a 
holographic renormalization procedure which we won't discuss in detail in this paper. 
Since the bulk fermions we consider don't couple to the dilaton and F fields directly, the 
results we obtain are still credible, in the absence of a full treatment of the holographic 
EMD theory. 



3. Holographic fermion with magnetic dipole coupling 

In order to explore the effects of magnetic dipole coupling on the spectral function of 
fermions, we start from the following action 

s f [$] = % J d d+2 x^ ^iv a v a — m — ipu)v + s hdy m (3.1) 

S bdy m = % j& d+l x^l^J^ + m_ (3.2) 



where S b dy is a boundary action to ensure a well defined variational principle for 
the total fermion action. And vj7 = vpr*, V a = (e a ) M £>^, with = d^-iqB^ + \uj^ ab r ab , 
T ab = |[r a ,r b ]. u^ab is the spin connection 1-form and H = ^T ab (e a )^(eb) u H^ u . T a are 
the d + 2 dimensional gamma matrices and (e a ) M are the vielbeins and M is the mass 
of the probe fermion. Furthermore, g e is the determinant of the induced metric on the 
constant r slice, r = e. ty± is defined by 

#± = -(i ± r r )# , rm ± = ±* ± (3.3) 

The Dirac equation derived from the action reads 

(T a V a -M- ipU)^ = (3.4) 
Taking a Fourier transformation 

*(r, x,) = {-gg rr )--*e- Mk ^^{r, fc„) , k^ = (-u, k) (3.5) 

where the prefactor was introduced to remove the spin connection in the equations 
of motion. Since the theory is rotational invariant, we can choose the momentum along 
the X\ direction. The Gamma matrices are chosen as follows 

r-c-f-y. M ! o\:.)' Mt-y ™ 

where a are Pauli matrices. We further set 
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Since the Dirac equation is first order, there exists some relation between ip + and 
if)-. Assuming if>+(r, k^) = — z£(r, kfj)if)-{r, k^), we can derive an elegant equation to 
extract correlators 

V^r£± + 2M£ d = ( v _ ± k^T^)e ± + {v+ =F k^/r^) (3.8) 
where £ + = iu-/u + , = id-/d + , C,± are the eigenvalues of the matrices £. And 

«± = V^iui + qB t ± P Vg^d r B t ) (3.9) 
The corresponding retarded functions can be readily obtained as follows [^6|] 

G (Ag = lim£(r,/g (3.10) 
At the event horizon, we impose in-falling boundary conditions 

£(r&,Ag=i, forw^O (3.11) 

We emphasize that the dimension of the fermionic operator O is A = (m + d)/2 
which leads to the unitarity bound was automatically satisfied with m > given by 



the null energy conditions [25|. The fermion mass decouples from the operator UV 
dimension, contributing only to the IR physics which is peculiar in the asymptotical 
geometries with hyperscaling violation. 

4. Numerical results and Emergence of the gap 

To extract the effects of bulk dipole coupling on the spectral function, we need to 
numerically solve the flow equation ( p.8|) with initial conditions ( |3.11|) . The spec- 



tral function is proportional to ImG(u,k), up to normalization. Due to the relation 
Gn(u, k) = G 2 2(u, —k), we will only consider G 2 2(u, k) and omit the subscript in the 
following. For convenience to perform numerical calculation, we set M — 0, \i — 1, q — 
2, z — 2, d — 3. The dipole interaction strength p and hyperscaling violation 9 (or n) 
remain to be free. 

First, we fix hyperscaling violation, considering n = case. From the plots above 
in figure [I], we find a sharp quasi-particle like peak at kp w 1.2044 for p — 0, indicating 
a Fermi surface. The dual liquid is of Fermi type with linear dispersion relation at the 
maximum height of the spectral function. For larger charge q, it has been investigated 
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Figure 1: The 3D and density plots of ImG{ui, k). In the plots above, p = 0, a sharp peak 
occurs at k ~ 1.2044. In the plots below, p = 4, a gap emerges around oj = 0. 



with great detail in |26] that more branches of Fermi surfaces will emerge and even 
a peculiar Fermi shell-like structure, which exactly contains many sharp and singular 
peaks in some narrow interval of the momentum space exists when q is large enough. 

In the plots below of figure 1, a gap emerges when the interaction strength was 
turned on p = 4 and there are two bands, located at positive frequency (we call it upper 
band) and negative frequency (called lower band) axis respectively. Interestingly, the 
upper band is much sharper than the lower one. For bigger momentum, the lower band 
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Figure 2: The 3D plots and density plots of ImG(uj,k) for p = 6 (plots above)and p = 8 
(plots below). 

is also as sharp as the upper band but becomes widened for relative small momentum, 
probably implying that the intensity of the spectral function tends to distribute at 
small momentum region in the negative frequency region. 

From plots in figure we can see that when p increases further, the gap becomes 
larger. The upper band always keeps sharp for all momentum, just translationally 
moving to the higher frequency region while the lower band is deformed much evidently 

■"■The black part in the 3D plots are purely numerical noise. 
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Figure 3: The plots of ImG{oj,k) for p = 0.01,0.2,0.4,1.1,1.2,2.0, from left to right and 
from top to down, respectively. 

by transfer of the spectral weight dominantly over the small momentum and negative 
high frequency region such that the intensity in the big momentum almost disappears. 

To explore further the distribution of the spectral function in detail when p in- 
creases, we show the plots of spectral function in figure [|. For very small p, the 
spectral function still has a sharp peak at u = 0, showing the main feature of a Fermi 
surface. As p increases, the intensity of the peak degrades and first transfers to nega- 
tive frequency region but soon returns to the positive frequency region when p increases 
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Figure 4: The gap width A as a function of p. 

further. Finally, around p cr u ~ 1.2 the original sharp peak at uj = disappears and two 
stable bands structure emerges in both positive and negative frequency regions. We 
observe that the dominant spectral density appears at the upper band. Furthermore, 
we show in figure |] that the width of the gap A becomes wider as interaction strength 
p increases. 

We now turn to vary hyperscaling violation 9 with fixed p. Without loss of gener- 
ality, we set p = 6. In figure |[ we show the 3D plots and density plots of the spectral 
function for n = 0.5 in the plots above and n — 1 in the plots below respectively. 
The gap width keeps almost independent of the hyperscaling violation while the bands 
structure are gradually smoothed out as n increases. We clearly find that the sharper 
upper band disappears much easier than the lower band. The similar phenomenon also 
happens to the lower band itself. As n is amplified, the relative sharper region of the 
lower band in high momentum disperses and becomes wider while the original rough 
region in the small momentum space becomes narrower and narrower, implying the 
effects of smoothing. We may argue that in the 9 — )■ d limit, any sharp peak of the 
spectral function will be completely smoothed out and the spectral density will transfer 
and redistribute to all frequency-momentum space homogeneously, with no explicit gap 
and bands structure again, probably indicating a critical phase. 

5. Conclusions 

In this paper, we have studied the novel features of fermions in the presence of bulk 
dipole coupling in the geometries with hyperscaling violation. For a finite hyperscaling 
violation 9 = d/2, we observe that when the dipole interaction strength p=0, a sharp 
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Figure 5: The 3D plots and density plots of ImG(uj, k) for n = 0.5, 1 in the plots above 
and below, respectively. 

quasi-particle like peak occurs near hp ~ 1.2044 at zero frequency, indicating a Fermi 
surface. As p increases, the intensity of the sharp peak degrades and first transfers 
to negative frequency region but soon redistributes and returns to positive frequency 
space. When p crosses a critical value p cr u, the Fermi sea disappears and a stable 
gap with two bands located on both positive and negative frequency axis respectively 
emerges for all momentum. It seems very interesting that these two bands behaves 
quite differently. The upper band on the positive frequency region keeps sharp for all 
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momentum while the lower band on the negative frequency region is only sharp for high 
momentum but disperses in the small momentum space. As p increases further, the gap 
becomes wider and the upper band is still sharp but the lower band redistributes the 
intensity of the spectral function to smaller momentum space by sacrificing the sharp 
region of the high momentum. We also find that the width of the gap increases with 
the increasing of p. 

When we fix p = 6 and turn on larger hyperscaling violation n = 0.5, 1 respectively, 
the gap almost doesn't change but the peaks and bands become smooth gradually. 
More interesting, the relative sharper upper band disappears first while the lower band 
also becomes narrower and narrower. Thus, the strength of the spectral density might 
distribute homogeneously to all frequency- momentum space in the 9 — >■ d limit. It is 
of certain interests to explore the postulated critical phase in this limit. We leave it in 
near future. 
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